Bonn-TH-2013-06 



SO{10) inspired GMSB 



o 

(N 

Oh 
< 

(N 
Oh 



> 

On 

O 
^* 

o 
m 

> 

X 
J3 



M. E. Krausfl and W. PorocQ 

Institut fiir Theoretische Physik und Astronomie, Umversitdt Wiirzburg 

Am Hubland, 97074 Wiirzburg, Germany 

F. StaulS 
Bethe Center for Theoretical Physics & Physikalisches Institut der Universitdt Bonn, 

53115 Bonn, Germany 

We consider a supersymmetric model motivated by a 50(10) GUT theory: the gauge sector 
near the SUSY scale consists of 5'L''(3)c x SU{2)l x U{1)r x U{1)b-l- We embed this model in 
minimal gauge mediation and incorporate neutrino data via an inverse seesaw mechanism. Also 
in this restricted model the additional D-terms can rise the light Higgs mass in a sizable way. 
Therefore, it is much easier to obtain rrih — 125 GeV without the need to push the SUSY spectrum 
to extremely large values as it happens in models with MSSM particle content only. We show 
that this model predicts a diphoton rate of the Higgs equal or smaller than the standard model 
expectation. We discuss briefly the collider phenomenology with a particular focus on the next to 
lightest supersymmetric particle where this model offers the sneutrino as an additional possiblity. 
Moreover we point out that also in this model variant supersymmetry can be discovered in Z' decays 
even in scenarios where the strongly interacting particles are too heavy to be produced at a sizable 
rate at LHC with 14 TeV. In addition, we show that lepton flavor violating observables constrain 
the size of the neutrino Yukawa couplings where in particular muon decays and fi — e conversion 
in heavy atoms are of particular importance. Once these constraints are fuUfilled the rates for r 
decays are predicted to be below the reach of near future experiments. 
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I. INTRODUCTION 

The LHC is rapidly extending our knowledge of the 
TeV scale. However, there is currently no hint of new 
physics beyond the standard model (SM) what leads to 
severe lower limits on the mass of new, especially col- 
ored, particles. One of the most popular model classes to 
extend the SM is supersymmetry (SUSY), in particular 
the minimal supersymmetric standard model (MSSM). 
As the MSSM itself has over 100 free parameters which 
are mainly part of the SUSY breaking sector, mechanisms 
of SUSY breaking have been studied which depend only 
on a few parameters and which predict distinct relations 
among the different soft terms. Those models trigger 
SUSY breaking in our visible sector by communicating 
with a hidden sector where SUSY gets broken at the 
first place. Popular mechanisms of transmitting SUSY 
breaking from the hidden to the visible sector work ei- 
ther via gravity like in supergravity [U [2| or via gauge 
interactions with so-called messenger fields like in gauge 
mediated supersymmetry breaking (GMSB) [5HTT] with 
six parameters. GMSB has the appealing feature that it 
is a flavor blind SUSY breaking. Hence, it solves auto- 
matically the flavor problem if the SUSY breaking scale 
is not too high. In addition, the gravitino is the light- 
est supersymmetric particle (LSP) and usually the dark 
matter candidate in this kind of models. 



However, both models are under big pressure because 
of the observation of a SM-like Higgs boson with mass of 
125 GeV [12l [13]. Even if this mass is below the abso- 
lute upper limit of about 132 GeV which can be reached 
for the light Higgs mass in the most general MSSM |14j . 
it is already very hard to explain it in the constrained 
MSSM (CMSSM) with five parameters ^Ml\ and de- 
mands large SUSY breaking masses and especially a large 
mass splitting in the stop sector. This mass splitting is 
caused by large trilinear couplings. However, in GMSB 
these terms are always small even if the ^i/Bf^ problem of 
the GMSB is solved HE]. That makes it even more diffi- 
cult to obtain a Higgs mass in the correct range TSj . This 
caused increasing interest in non-minimal GMSB models 
which involve also superpotential interactions between 
the matter and messenger sector to create large trilinear 
terms 
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Another possibility to reduce this tension between the 
Higgs mass and the simplest constrained models is to ex- 
tend the Higgs sector of the MSSM. The smallest possible 
extension is to add a gauge singlet like in the next-to- 
minimal supersymmetric standard model (NMSSM) (see 
p8l [29] and references therein) . The singlet and the cor- 
responding superpotential coupling to the Higgs doublets 
can significantly lift the upper limit on the light Higgs 
mass of rrih < Mz at tree-level in the MSSM by new F- 
term contributions [30 . It has been shown that even in 
the constrained NMSSM a Higgs mass of 125 GeV can be 
explained [311 132] . If one drops in addition the assump- 
tion of a Z3 symmetry and consideres instead the gener- 
alized NMSSM these masses are obtained with even less 
fine-tuning [33,. The same feature can also be observed 



in models with Dirac instead of Majorana gauginos which 
usually come with an extended Higgs sector [34j. There 
are also some hints that the branching ratio of the ob- 
served particles do not completely agree with the SM 
expectations. Especially, the diphoton rate seems to be 
enhanced what is usually hard to explain in the context 
of the (constrained) MSSM |3Fl[35]. In contrast, this en- 
hancement can much easier be obtained in the NMSSM 
or its generalized version |3TJ ISTtiiO] . 

A second possibility is to consider models with ex- 
tended gauge structures which arise naturally in the con- 
text of embedding the SM gauge group in a larger group 
such as 50(10) or Eq, see e.g. [HM5] . In those models 
the upper bound on the lightest Higgs boson is also re- 
laxed, this time due to additional i'-term contributions 
[47H52] . Furthermore, they often provide the possibil- 
ity to explain neutrino data: either via seesaw type I-III 
which involve heavy states J53H56J or via inverse or linear 
seesaw with additional matter fields at the SUSY scale 
[57H59] . Another interesting feature is that such mod- 
els can also have a new gauge boson {Z') with a mass 
in the TeV range (SOI IM]- Hence, intensive searches for 
Z' bosons have been performed and bounds on its mass 
have been set [621 - I65] . For reviews on various Z' models 
see e.g. [551157] . 

In addition, U{1) extensions of the SM provide another 
peculiar feature which can have very interesting effects, 
namely gauge kinetic mixing [55H70) . While gauge ki- 
netic mixing is often ignored in phenomenological studies 
it has been shown recently in several works that it can 
have a significant effect on Z' phenomenology [7TH73], 
but also on the Higgs mass [75 and dark matter proper- 
ties [HI [76] . 

In this work we will assume minimal GMSB inspired 
by SO{10): the GUT group gets broken to SU{3)c x 
SU{2)l X U{1)r X U{1)b-l very close to the GUT scale 
but well above the scale of the messenger fields which 
trigger GMSB. In contrast to previous studies, we as- 
sume however, that the breaking down to the SM gauge 
groups takes place at the TeV scale, i.e. well below the 
lowest messenger scale. Hence, we study a messenger 
sector charged under U{l)iiX U{1)b-l what will change 
our boundary conditions. Furthermore, we gain an en- 
hancement of the mass of the light Higgs boson. In ad- 
dition, we follow the setup of [77] where this model has 
been studied in gravity mediation and assume additional 
gauge singlets present at the SUSY scale to incorporate 
neutrino masses and mixing via inverse seesaw. 

The remainder of this paper is organized as follows: 
in Sec. [IT] we briefly summarize the main features of the 
model and its particle content. In Sec. [TIT] we discuss 
the results for Higgs physics under considerations, check 
the bounds coming from LFV observables on the model 
parameters and comment on the expected collider phe- 



II. ASPECTS OF THE MODEL 

A. Particle content and superpotential 

In this section we discuss briefly the particle content 
and the superpotential of the model under consideration. 
For a detailed discussion of the particle spectrum we refer 
to [521 EZ] • The superpotential is given by 



+ Yl^i'-LjHu + Y'g'v^S.jXR - fJ.RXRXR + fJ^'s^zSj , 

(1) 

where the upper line corresponds to the standard MSSM 
superpotential and the lower line contains the new sector 
as well as the ingredients for the inverse seesaw mecha- 
nism: Ys and Y^ being the neutrino Yukawa couplings 
and fis the mass term for the singlet field S which is re- 
sponsible for the mass of the light neutrinos. In Sec. |IID| 
we will see that this model with an inverse seesaw mecha- 
nism for neutrinos is much easier to implement in GMSB 
than the corresponding model with type I seesaw where 
the Si fields are absent. 



The scalar fields XR and XR break U{1)r x U{1)b-l 
down to U{1)y while conserving i?-parity. However, to 
ensure this an extra Z^ matter parity has to be intro- 
duced, such that not all terms allowed by gauge invari- 
ance are present in eq. ([I]). Interestingly, the particle 
content of this model is in agreement with gauge cou- 
pling unification even if the breaking scale of SU{3)c x 
5C/(2)l X U{1)r X U{1)b-l ^ SU{3)c x SU{2)l x U{1)y 
is close to the breaking scale down to SU{3)c x U{l)em- 
Therefore, we will always assume a one step breaking 
SU{3),xSU{2)lxU{1)rxU{1)b-l ^ SU{3),xU{l),^ 
in the following. However, to facilitate the comparison 
with the MSSM we will work in a different basis for the 
U{1) sector: we will take t/(l)y x [/(l)^ as the orthog- 
onal basis instead of U{1)r x U{1)b-l- We give the 
corresponding U{1) quantum numbers for both bases in 
Tab. [I] 

The soft SUSY breaking terms are 



1 



nomenology. We conclude in Sec. IV 



Vsoft = mij(j)*(l)j + ( -MabXa>^b + Bf^HuHd + Bf^^XRXR 

+ Tl^HaefL, -f T^^H^i^^L, + Tl'xR^^SJ + h.c. 

(2) 

with the generation indices i and j . We have introduced 
here (t)i for all scalar particles and Xa for the different 
gauginos. Note that, because of the two Abelian gauge 
groups present in the model and the consequential gauge 
kinetic mixing discussed in the next subsection, also the 
mixed soft gaugino term Myx^yX^ is present j78] . 
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TABLE I. Chiral superfields and their quantum numbers with respect to S'f/(3)c x SU{2)l x U{1)r x U{1)b-l- We also give 



the quantum numbers in the basis SU{Z)c x SU{2)l x U{1)y x (7(1);^, the relations between both bases is defined in Sec. II B 



B. Gauge kinetic mixing 

Even if t/(l)fl and U{1)b-l can be embedded orthog- 
onal in S'O(IO) at a given scale a kinetic mixing term of 
the form 



Lr 



= -xF 



B-L 



fii^pl 



(3) 



can occur. The reason is that the Higgs fields we assume 
to be present at the SUSY scale do not form a complete 
representation of 5*0(10). Hence kinetic mixing will be 
introduced by RGE evolution. This can be seen by the 
off-diagonal elements of the anomalous dimension matrix 
which in the basis {U{1)r,U{1)b-l) at 1-loop is given 
by 



(4) 



Here A'' = diag(l, \/3/2) contains the GUT normaliza- 
tion. In order to correctly account for gauge kinetic mix- 
ing effects we follow the approach given in |75] and shift 
the term to a covariant derivative of the form 




D^'=d^- iQiG. 
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where 



G = 



9r 9rbl \ 
^9blb. 9bl ) 



(5) 



(6) 



A^ = {A^^A^_j^)'^ and Q is a vector containing the 
[/(I) charges of the field under consideration. We assume 
the breaking into U{\)ii x U{1)b-l to take place at the 
GUT scale Mqjjt and demand 9!iBL — 9blr = at 
Mgut as initial condition. In addition, we have the free- 
dom to go into a particular basis by rotating the gauge 
bosons of the Abelian groups. As already mentioned, 
for an easier comparison with the usual GMSB we take 
the basis U{1)y x U{l)y^ where the first factor is the 



usual hypercharge and the second one is the orthogonal 
one within 50(10). The gauge couplings and charges of 
U{1)r x U{1)b^l and [/(l)y x U{1)^ are (without GUT 
normalization) related via 



A^ ^ A'f" 



G^G' = 




Qb-l + Qr 

hB-L-QR 



(7) 



with 
9Y = 



9bl9r — 9blr9rbl 
\/{9blr - 9R.y 



9BL - 9rblY 



5x = 5 \/{9BLR - 9rY + {9BL - 9rblY ■ 

9yx = 
2(5|l 



9%Lr) + 9BLR9R + 9BL9RBL ~ ^9^ + 9rbl) 



^VidBLR - 9rY + {9BL - gRBbY 



C. GMSB boundary conditions 



(8) 



In GMSB models it is assumed that supersymmetry 
breaking is generated by one or more superfields X^ hv- 
ing in a "secluded" sector. We assume for simplicity that 
only one field X is present which is coupled to a set of 
messenger superfields ^^ via 



Wgm = KX $,$, 



(9) 



Furthermore, it is assumed that the scalar and auxiliary 
components of X receive a vev 



{X) = M + e^F , 



(10) 



and that it couples universal to $i implying that one 
can set Ai = 1. The supersymmetry breaking due to 
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which can be estimated to be 



TABLE II. Quantum numbers of the messenger fields in the 
respective bases. 



the F-term vev is communicated to the visible sector 
via the gauge interactions of the $i. Since we are inter- 
ested in minimal gauge mediation without spoiling gauge 
coupling unification we assume that the messenger fields 
form a complete 5*0(10) multiplet, e.g. a 10-plet. This 
results in two SU{2)l doublets and two SU{3)c triplets 
below the 5*0(10) scale with suitable charges under the 
Abelian gauge groups which are listed in Tab. [llj 

The SUSY breaking gaugino and scalar masses are gen- 
erated via 1- and 2-loop diagrams, respectively 79 - 8^ . 
Neglecting gauge kinetic mixing the boundary conditions 
for the SUSY breaking masses are given by [5T] 
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(11) 



(12) 



with A = F/M, Xi — |A/M|, g{x) and f{x) are approx- 
imately 1 for a; < 0.2. ga denotes the coupling of gauge 
group a and i runs over the messenger fields. na{i) is the 
Dynkin index of the messenger with respect to the gauge 
group i. We use a normalization where n^ = 1 for the 
10 of 50(10). Ca{k) is the quadratic Casimir invariant 
of the scalar field k. 

For a proper treatment of gauge kinetic mixing we use 
the substitution rules for Abelian groups given in [75]. 
The resulting soft masses for the gauginos and scalars at 
the messenger scale read 
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(14) 
^-(16^^'( E CA{k)g\Y.fi^^)^A{^) 

^ ^ A^Abelian i 

+ J2 f{xi){QlNGG'^NQ,f) . (15) 



The trilinear soft SUSY breaking parameters are, as 
usual in minimal GMSB, essentially zero at the scale of 
gauge mediation. The singlet 5* is a special case because 
it is a gauge singlet and, thus, would have a zero mass 
at this level. However, it gets a mass at the 3-loop level 



y2 



(16) 



Obviously this mass squared parameter is suppressed by 
an additional loop factor and RGB effects usually drive 
it to negative values at the elcctroweak scale. However, 
as can be seen in Sec. |HIC3| this is compensated by 
an _F-term proportional to AI^, yielding a positive mass 
squared for the corresponding mass eigenstates. 

For completeness we note that one can explain neutrino 
by adjusting fig and taking Y^, as well as Yg diagonal. 
fis is a small parameter which does not affect the collider 



phenomenology. However, we will discuss in Sec. HI E the 
effect of non-diagonal entries in Y^, and Yg in the range 
compatible with neutrino data on rare lepton decays. 



D. Tadpole equations 

We decompose the neutral scalar fields responsible for 
gauge symmetry breaking as usual: 
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(17) 



We use the minimization conditions to determine the pa- 
rameters l/zp, l^iflP, Bfj, and -B^,^: 
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where gL is the SU{2)l gauge coupHng, t^, c^, Sx 
tana;, cos a:, sin a; whereas tan/3 = ^^, tan/Sn 



vl and v\ 



corresponding terms can be generated by the Giudice- 
Masiero mechanism |83| and are thus free parameters in 
our context. 

The latter of these equations is of particular interest 
as it is responsible for one of the major limitations to 
the model. The 1-loop /3 functions for the soft-breaking 
masses read in the limit of vanishing kinetic mixing 



Note, that the singlet field 5'. Technically this amounts in replacing the 
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(24) 



which is zero at the messenger scale and which stays zero 
if only 1-loop RGEs are used. One can see that the main 
differences in the running are stemming from terms which 
are proportional to the trilinear soft-breaking couplings 
or the soft-breaking masses. Since we will consider the 
minimal GMSB where non-vanishing trilinear couplings 
are only generated via RGE evolution and the breaking 
takes place well below the GUT scale, the splitting be- 
tween the soft parameters m?, and m% will in general 
be smaller in comparison to a scenario with gravity me- 
diation. Because of eq. (21) this immediately constrains 
tan/3fl; to be larger than but close to one. The terms 
proportional to the veva squared then only give nega- 
tive contributions to l/x^p, i.e. there is an upper limit 
on \vr\ depending on tan/3/j to find a solution to the 
tadpole equations. In Figure^ we show j/i^l in the vr- 
tan/3/j plane where one can see the correlation between 
the two parameters. Note, that in the upper white area 
one cannot achieve the correct gauge symmetry breaking 
whereas in the lower white area one encounters tachyonic 
states. 

We have also considered the case where neutrino 
masses are generated by a seesaw type I mechanism simi- 
lar to [74 , in which case there is no need to introduce the 



terms Y^' ^f S, XR + Ms' S^ S, in eq. (jlf by rf ^>f Xfl^| 
where x'r has twice the U{1) charges of xr- Perform- 
ing the same chain of calculations one finds that there 
are hardly points with broken C/(l)x as the larger gauge 
contributions in the RGE evolutions prevents tti^' to be- 
come sufficiently small. 
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FIG. 1. Allowed parameter space in the vr — tan/3_R plane. 
The plotted values correspond to \fJ,R\ which is calculated us- 
ing the tadpole equations. The free parameters have been 
set to n = 1, A = 5 • 10^ GeV, M = lO" GeV, tan/? = 
30, sign(pfl) = -, diag{Ys) = (0.7,0.6,0.6), F," = 0.01. 



E. Higgs sector 

In GMSB models with MSSM particle content one 
needs a SUSY spectrum in the multi-TeV region to ac- 
commodate a Higgs mass of 125 GeV (see, e.g., P^ [Ml 
[55]). The reason is that the trilinear soft SUSY breaking 
couplings are zero at the messenger scale in the minimal 
model and thus the loop corrections to the Higgs boson 
masses get reduced. In our model the additional U{1) 
factor gives already a sizable Z3-term contribution to the 
tree-level part of the Higgs mass and, thus, the need for 
large loop corrections gets reduced. 

On tree level, the scalar Higgs mass matrix in the basis 
{(7d,(Ju, o-R, (Tr) is given by 



mlo = 








f \~9h^cl -h m\sl 


-^-f(gy^+4m\) 


l~glvvRCfiCf}^ 


-lglvvRCfiS(i^ 




-'-f{~9lv^ + ^m\) 


i9h^4 + "^'acj 


-IglvvRSpcp^ 


yivvRspsp^ 




l~glvvRcpcp^ 


-l~glvvRspCfi^ 


Wx^^Ir+'^Ir^r 


-^(255X + 16miJ 




\ -UlvVRCfiSfi^ 


IglvVRSfsSf!^ 


-^(25ff2 4 + 16miJ 


25 2 2 2 ,2 2 

iqBx^r'^Pr ^^ArCp^ 



(25) 



where g^ ^ gl +g^ + (g^ - gy^f, g^ = g^{g^ - gvx) 
and Sx^ Cx = sinx, cosx. The parameters mA and 
toah are the tree-level masses of the pseudoscalar Higgs 
bosons, which are given by m\ = B^/{spCp) and m\ — 



Btj,f./{sj3^ci3j^). Already at tree level this leads to non- 
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FIG. 2. Tree-level dependence of the lightest Higgs masses 
(above) as well as the admixture of the SU{2)l doublet Hig- 
gses 7?. I,. = |t/ii|^ + |(7i2p (below) on tan/3_R with the param- 
eter choice of Fig. [l] and vr = 7 TeV. The horizontal small 
dashed (red) line shows the Z mass. 



F. Dark matter 

As already mentioned, the gravitino is the LSP in 
GMSB models and all SUSY particles decay into it in 
a cosmologically short time [55M88J . The abundance of 
thermally produced gravitinos is under assumptions con- 
sistent with the standard thermal evolution of the early 
universe given by 



^3/2^ — 



2 _ "I3/2 100 



keV 



(27) 



Here g* is the effective number of degrees of freedom at 
the time of gravitino decoupling. For a mass of O(IOO) eV 
the gravitino would form warm dark matter and would 
have the correct abundance to explain the observed dark 
matter relic density in the universe. However, there are 
stringent constraints on the contribution of warm dark 
matter from observations of the Lyman-a forest |89) . 
These bounds rule out pure warm dark matter scenarios 
with particle masses below 8 keV for non resonantly pro- 
duced dark matter [SD] . If one takes this lower limit into 
account one sees that gravitinos which have once been in 
thermal equilibrium would overdose the universe. This 
is known as the cosmological gravitino problem. There 
have been some proposals in literature to circumvent this 
problem by, for instance, additional entropy production 
after the freeze out of the gravitino [9T1I93J . However, 
it turned out that entropy production from messenger 
decays hardly works [M]. Hence, one has to assume ei- 
ther other mechanisms like saxion decays [5S] or decays 
of moduli fields [96] • Also if the gravitino mass is in the 
MeV range they might never have been in thermal equi- 
librium if the reheating temperature is sufficiently low 
[57] . Because of these very model dependent issues we 
don't address the question of the gravitino relic density 
in the following. 



III. NUMERICAL RESULTS 



negligible contributions to the doublet Higgs mass as 
ig|i;2 ~ A/| -I- Iv^ig^ - gy^)^ > M|. As typical val- 
ues we find g^ — gy-^ ~ 0.27. This immediately gives an 
upper bound on the tree-level Higgs mass: 



mh,tree < Mz + T (Sx ^ 9Yx) 



(26) 



In Figure [2] we show the dependence of the two lightest 
Higgs states on tan/J^ at tree level. Even in this very 
restricted model the tree-level mass can easily reach 100 
GeV while at the same time requiring that this state is 
mainly an SU{2)l doublet Higgs boson. Even though 
the details are changed by loop corrections, see e.g. [5^ . 
this figure also shows that tan/3fl has to be close to 1 
to obtain this desired feature. In the numerical part we 
include the complete 1-loop correction to eq. ( 25 1 and the 



A. Implementation in SARAH and SPheno 

We used the implementation of the model in SARAH 
PHTUT] and SPheno [TUH HUS] presented in [77] and ex- 
tended it by the GMSB boundary conditions: here we 
allow for up to four messenger 10-plets with degener- 
ated masses. At the messenger scale we implemented the 
GMSB boundary conditions for the soft masses using the 
expressions eqs. (13)-(15). The link between SARAH and 



dominant 2-loop corrections to the MSSM sub-block. 



SPheno allows for a precise mass spectrum calculation 
based on full 2-loop RGE running and the 1-loop cor- 
rections to all masses. In addition, the known 2-loop 
corrections to the Higgs masses in the MSSM are linked 
|104H106] . For more details of the mass spectrum cal- 
culation as well as the inclusion of SUSY thresholds we 
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A [GeV] 
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- 




+ 
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Wfl [TeV] 
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0.01 




0.01 


0.01 


0.01 


djap(ys) 


(0.65,0.65,0.1) 


(0.65,0.65,0.3) 


(0.65,0.65,0.65) 


(0.6,0.6,0.6) 


(0.77,0.73,0.45) 


(0.7,0.6,0.6) 


m^, [GeV] 


70 


92 


125 


70 


108 


98 


^i.'.i 


0.006 


0.018 


0.961 


0.003 


0.094 


0.006 


m^, [GeV] 


126 


127 


156 


124 


124 


124 


^i... 


0.994 


0.982 


0.039 


0.997 


0.906 


0.995 


M^, [TeV] 


2.53 


2.7 


2.41 


4.32 


m.,,, [GeV] 


357 


1070 


2306 


2277 


1542 


3633 


m.,,, [GeV] 


2309 


2308 


2306 


2278 


2497 


3633 


'""ft. 3 [GeV] 


2309 


2308 


2306 


2278 


2633 


4238 


msi [GeV] 


334 


909 


1715 


1728 


1207 


1863 


mi>2 [GeV] 


1072 


1546 


1715 


1757 


1482 


1879 


ms3 [GeV] 


2090 


2048 


1715 


1759 


1514 


1879 


mi=i [GeV] 


906 


906 


905 


867 


764 


1007 


m^^ [GeV] 


1166 


1166 


1165 


976 


877 


1061 


me-« [GeV] 


1167 


1166 


1166 


976 


877 


1061 


m^o [GeV] 


505 


766 


1156 


575 


453 


589 


m^o [GeV] 


1157 


1157 


1353 


610 


825 


1043 


m^± [GeV] 


2216 


2216 


2217 


1113 


883 


1142 


m,± [GeV] 


2591 


2590 


2588 


1956 


1600 


2015 


mg [GeV] 


5460 


5459 


5456 


3018 


2423 


3076 


m,, [GeV] 


4209 


4209 


4206 


2993 


2231 


2941 



TABLE III. Input parameters and mass spectrum of different representative parameter points. 



refer to [77|. In addition, the SPheno version created by 
SARAH includes also routines for a full 1-loop calculation 
of the LFV observables k — ?► Ijj, U -^ ilj, (x — e conver- 
sion in atoms, flavor violating r decays to a lepton and 
meson as well as Bg -^ l^^lJ-~ |107j . 

For further discussions we choose six benchmark sce- 
narios BLRI - BLRVI which provide distinct features. 



These benchmark points are given in Tab. Ill and will be 
discussed in the next following subsections. 



B. Higgs physics 

We performed a scan over the free parameter space in 
order to numerically check how well Higgs data can be 
accommodated for in our GMSB framework. The pa- 
rameter variations can be found in Tab. |IV[ In Fig. [3] 
we show the masses of the doublet-like Higgs versus the 
mass of the lightest stop. 



Parameter 


varied range 


# Messenger multiplets n 


1 ... 4 


Messenger scale M 


(10^ ... 10^2) Qgv 


A^ F/M 


^(10^ ... 10«) GeV 


tan/? 


1.5 ... 40 


tan^H 


1 ... 1.15 


sign(^H) 


±1 


Vr 


(6.5 ... 10) TeV 


y^' 


0.01 ... 0.8 


y:^ 


10-^ ... 0.5 



TABLE IV. Parameter ranges of the scan. The sign of /i has 
aways been taken positive. 



As expected from the discussion in Sec. |IIE[ points where 
/i2 is the doublet- like Higgs are of particular interest since 
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FIG. 3. Mass of the doublet-like Higgs versus the mass of the lightest stop for n — 1 and the other parameters as in Tab. |IV| 
Only points with Rh^-y^ > 0.5 (left) and 0.9 (right) were included. The blue dots represent points where the lightest eigenstate 
is doublet-like, green dots where it is the second-lightest Higgs. 



they allow for higher values of the Higgs mass at a fixed 
mf. Due to the tree level contributions from the new 
sector, we can achieve the observed Higgs mass even for 
stop masses of about 2 TeV while a doublet-like hi re- 
quires mf > 3 TeV. Admittedly, this is quite a high scale 
in terms of naturalness in SUSY. However, compared to 
the lower limit of rrij: > 5 TeV in usual GMSB scenarios 
(see, e.g., j$4)) this is significantly lower. Such a heavy 
stop will be difficult to study at the LHC and will poten- 
tially require a center of mass (c.m.s.) energy larger than 
14 TeV. However, here an e+e~ collider like CLIC with 
up to 5 TeV c.m.s. energy might be an ideal machine to 
discover and study such a heavy stop, see e.g. |108| I109J 
and refs. therein. 

A way to allow for a lighter SUSY spectrum in GMSB 
scenarios apart from the mixing with the extended Higgs 
sector is going up to higher messenger scales, thus allow- 
ing a longer RGE running and hence larger induced T- 
parameters as demonstrated in Fig. |4J At M ~ 10"'^° GeV 
and 10^^ GeV there is a level crossing between the two 
light states for the points BLRHI and BLRV, respec- 
tively, which is the reason for the observed increase of 

An interesting observable is the rate ft, — >■ 77 as there 
are some hints for an enhancement above SM expecta- 
tions |110H112] . We define the ratio Rh^^y^ by 



Rh- 



m 



[ajpp -^ h) X BR{h -^ J'^)]blb. 
[a{pp -^ h) X BR{h -^ 77)]sM 



(28) 



The cross sections for the main production channels, 
gluon fusion and vector boson fusion, are essentially the 
SM-production cross section re-weighted by the (effec- 
tive) couplings of the Higgs boson cf^x^ normalized to 
the SM expectations c^xx'- 
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FIG. 4. r„ and mass of the doublet-like Higgs versus M 
for BLRIII (black dashed line), BLRIV (black full line) and 
BLRV (black dotted line) (but for tan/3_R = 1.03). The hght 
green lines correspond to the Xfi-like Higgs state of the cor- 
responding parameter point. 
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FIG. 5. D-term contribution to the mass entries of the 
i?-sleptons for tan/3 = 10, Mz' = 1.5,3.0 TeV and fix- 
ing the gauge couphngs by the requirement of gauge cou- 
pling unification: g^^ — 0.55, gy = 0.36 respectively 



^fixs-x. = 0.57, gn = 0.45. 
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The full (dashed) lines corre- 
L {U{1)y X U{1)b-l) scenario. 



Q^ 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 






600 



800 



1000 



1200 



1400 



[GeV] 



FIG. 6. Decay rate of h — > 77 as a function of the stau 
mass using four 10-plets. Only points with 123 GeV < rrih < 
128 GeV are included. Color coding of the parameter points 
as in Fig. [3] 



C. Next to lightest SUSY particle 
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(29) 



The main contribution to Higgs production comes from 
gluon fusion. The effective Higgs couphng to two gluons 
is completely determined in the SM by the top and W 
loop. 

In supersymmetric models an enhancement can be 
achieved via a hght stau. In models with extended gauge 
structures such a light stau and thus an enhancement of 
the 77 rate can be obtained even in scenarios with large 
soft SUSY breaking parameters |113j as there are large 
negative contributions due to the D-terms of the extra 
U{1) to the stau mass. However, in the model consid- 
ered here this does not work for two reasons: the large 
stop mass required to obtain the correct Higgs mass im- 
plies a lower limit on A. Secondly the Z3-term itself is 
smaller in our model compared to the one of ref. [113] 
taking the same Z' mass and ratio of additional vevs as 
demonstrated in Fig. [5J In this restricted model the tad- 
pole equations imply that a larger vr requires a smaller 
tan/3jj and, thus, the I?-terms cannot be enhanced to the 
required level. In Fig. 6] we show Rh^jj as a function of 
the stau mass demanding that the mass of the doublet- 
like Higgs to be in the range 123 GeV < rn/i < 128 GeV. 
This implies a lower limit on TTif ^ > 500 GeV which is too 
large to get a sizable contribution to ft. — >■ 77 and thus 
we find Rh^^j < 1 in this model. Hence, this model will 
be excluded if Rh^-yy > 1 is established by ATLAS and 
CMS at a significant level. However, the most recent 
results of CMS point exactly in this direction that the 
diphoton rate is in good agreement with SM expectations 
[114, . 



As usual in GMSB models the gravitino G is the light- 
est supersymmetric particle (LSP) and its mass is given 

by ISO 



W3/2 



\/3; 



mpi 



(30) 



with the reduced Planck mass mpi . The gravitino mass is 
usually in the MeV range or above due to large messenger 
scales required to obtain the correct symmetry breaking. 
In our model are three possibilities for the next-to light- 
est supersymmetric particle (NLSP): the two usual can- 
didates which are the lightest neutralino and the lightest 
slepton, which is usually a stau. In addition our model 
contains the lightest sneutrino as third candidate. In gen- 
eral we can state that the lightest neutralino will be the 
NLSP for low messenger multiplicities, n < 2, and little 
hierarchy in the diagonal entries of Yg. For larger n, Xi 
can only be lighter than the lightest slepton if the left- 
right-splitting of the stau is small {i.e. for low tan/3 val- 
ues) or if |/Xii;| is small. We discuss the different character 
of a neutralino NLSP in the next subsection. Otherwise, 
i.e. for large n and non-hierarchical Ys, the stau is the 
NLSP. A sneutrino can be the NLSP for all n if there 
is a large hierarchy in the Ys entries: the scalar singlet 
field corresponding to the smallest Y5 entry gets light. In 
the following we present the corresponding mass matri- 
ces and discuss briefly the main differences compared to 
the phenomenology of the usual minimal GMSB model 
using the parameter points in Tab. Hi] As the lifetime of 
the NLSP is proportional to F^ [5D] we find that in most 
of the available parameter space the NLSP is so long- 
lived that it will leave a typical collider detector before 
decaying. However, its lifetime is in general still below 
the bounds set by big bang nucleosynthesis. 
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1. Neutralinos 

This model contains seven neutralinos which are beside 
the usual MSSM gauginos and higgsinos the extra U{1) 



gaugino X^ and the two i?-higgsinos XR and XR- In the 
basis (Ay, Av^^, h^, h^, A^, XRi Xr) the mass matrix reads 



Mp 





gyVd 

1 

gyfu 
2 

2 







Af2 


ffi-fd 


2 
2 


QLVd 

2 
2 







(9x-9Vx)"d 


2 





2 
2 





2 





V 



2 



(3x-gyx)'"d 

2 

(3x-gyx)"" 

2 

, 4 
5gx"xa 









sgxi-x 



4 











Sgx^xfl 

4 
-fJ-R 




(31) 
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FIG. 7. Neutralino masses as a function of hr for the point 
BLRIV {Ml ~ 575 GeV) and 1.02 < tan/3fl < 1.033 to satisfy 
the tadpole equation (211. 



For a first understanding it is useful to neglect the mix- 
ing between the MSSM states and the additional ones. 
In this case one gets M^, ~ j^g^vj^ and in the limit 
tan/3/j — )■ 1 one finds for the eigenvalues of the three 
additional neutralino states 



fJ-R, 




^J■R± 



-A/2 
4 ^ 



Aq 



2M^^iR + ^2, 



(32) 



In most of the parameter space one finds |/ifl,|, M^ <C Mz' 
and thus one has one state with mass \^.r\ and two states 
with masses close to Mz' which can even form a quasi- 
Dirac state. For the MSSM-like states the lightest one 
is always bino-like and, thus, we find depending on the 
ratio \iir\/Mi that the lightest neutralino is either bino- 
like or a nearly maximal mixed XR ~ XR. state. This is 
exemplified in Fig. [7|for the point BLRIV with a slight 
adjustment of tan fi^ to satisfy the tadpole equation ( 21 1 . 




Here, the NLSP nature changes from Xi^-like to bino-like 



200 400 600 800 1000 1200 1400 

m [GeV] 

FIG. 8. Branching ratios of Xi (above) and of X2 (below) as 
a function of ^ir for the same parameter choice as in Fig. [T] 



at about ^xr ~ 575 GeV. 

The lightest neutralino will decay dominantly into a 
Xfl-like Higgs state and a gravitino G if it is mainly a xb_ 
higgsino whereas the MSSM-like bino state decays dom- 
inantly into 7G and ZG as depicted in Fig. [sj However, 
as mentioned above, the neutralinos are rather long-lived 
and, thus, at the LHC they will decay in general outside 
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the detectors. Hence, new techniques would be neces- 
sary to observe these states. For \fiii\ < Mi we find that 
h^^ can be produced in the decays of X2 ^* ^ sizable 
rate. Therefore SUSY cascade decays offer the possibil- 
ity to study this particle which can hardly be produced 
directly or in Higgs decays. 

The large values of A imply that the squarks and the 
gluino are usually in the niulti-TeV range implying that 
one will need the high luminosity option of LHC to study 
these particles in detail. It turns out that the two lightest 
states are g and ii. Depending on mg—mi_^ the gluino de- 
cays either dominantly into third generation quarks and 
neutralinos/charginos or into tti. In both cases the final 
states will contain 5-jets and W^-bosons. Depending on 
the nature of the two lightest neutralinos also a Higgs 
boson can be in the final state as discussed above. More- 
over, also the additional sneutrinos can appear in the 
cascade decays but distinguishing them from the usual 
MSSM sneutrinos will be rather difficult. 



2. Charged sleptons 



As explained above, stau NLSPs can be obtained for 
n > 3 and large values of tan /3 as left-right mixing can 
compensate the additional D-ierTus. As can be seen in 
Fig. |9] the three diagonal entries of I5 have to be of 
roughly the same size. 



As indicated in Tab. HI the gluino is usually very heavy 
and it turns out that it is the heaviest strongly inter- 
acting particle. This implies that one will need a very 
high luminosity to discover this particle. In general it 
decays into all squarks which in turn decay further into 
the MSSM-like neutralinos and charginos. Note, that 
both light Higgs states, the doublet-like one as well as 
the /i^jj-like one, can be produced in these decays. Fi- 
nally the lightest neutralino will decay into rfi and f 
will in general decay outside the detector. Thus, a typ- 
ical event will consist of several jets and leptons plus a 
charged track from a (at the detector-level) stable parti- 
cle. The phenomenology of long-lived staus has already 
been studied comprehensively in the literature, see, e.g., 
[1151 I116J , and bounds of ruf > 300 GeV have been set 
by the LHC collaborations [11 7j in MSSM scenarios. 



The mass matrix of the sleptons reads in the basis 
{eL,eR) 



(33) 



^{Tecp - ^i*Y,Sfi) ml + \v^clY,Y^ + DrI) ' 

(34) 

which has the same structure as in the MSSM but for the 
concrete form of the D-ievras 

Dl =^ (2( - igl + 9x9Yx + 2(.9y - s! + gl^)yc2p 
~ ^9xi'^9x + '^9yx)vrC2p^j , (35) 

Db =^ (2(.9x + ^9x9Yx - <9y + 9Yx)y'^2p 

+ ^9xi9x + 4gi'x)"flC2/3„ ) ■ (36) 



3. Sneutrinos 

As it is well-known, in inverse seesaw scenarios, the pa- 
rameter /i5 has to be small to explain correctly neutrino 
datej^ For completeness we note, that the inverse seesaw 
mechanism yields in this model three very light Majo- 
rana states which can explain the observed neutrino data 
where the six heavier neutrinos are pairwise degenerate 
forming three quasi-Dirac states [77] . We will denote the 
former by v and the later by vu ■ The i^-terms induced by 
fis as well as the corresponding soft SUSY-breaking term 
Bfj,^ SS induce a splitting of the complex sneutrino fields 
into their scalar and pseudoscalar components. However, 
in practice this mass splitting is tiny and, thus, we can 
safely neglect it in the following discussion. In the limit 
fiSj ^ps ~^ the sneutrino mass matrix reads in the basis 
(i>, i>S5) 



M? = 



( 2 



2 '^ V '^v 



D'l^ 



■^{TvSp - fi*Y^cp) mla 
V ^vvRY^Y^spSfi^ 



%/2 






(37) 



D'nl 



^vvRY^YsS/sSf)^ \ 
lAiYsCfii, 



%Tss,, 



'-^Y^Ys 



I 



^ For a discussion and the corresponding mass matrix see for ex- 
ample in | 77| and references therein. 
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FIG. 9. Masses of the lightest SUSY particles (left) and the sneutrino composition (right) as a function of Yg^ for the parameters 
specifying the points BLRI, BLRII and BLRIII. In the right plot the full red (dashed blue) line gives the 5* (D^) contribution 
to the nature of the lighest sneutrino. 



with 



D'l ^^[^i ~ 3ffx + 9x9Yx + "^{gl +9y+ 9Yx))'"^'^^P - ^9xi^9x + '^9Yx)'"RC2Pn 

D'r =-^{^(9x - 9Yxyc2i3 + 5gxw|fC2fe 



(38) 
(39) 



Obviously, the masses of the sneutrinos depend 
strongly on Ys- In particular S-S entries are dominated 

2 2 

by "^'^-" Yg Ys because to| is rather small as discussed in 



Sec. II C Therefore, even mi < does not automatically 



imply spontaneous i?-parity breaking. The iy'^-S mixing 
entry can be of the same size as the corresponding diag- 
nal entries for sufficiently large \fJ-R\. The entries which 
mix these states with i>i are much smaller and can be 
neglected for the moment. As tan/^/j is close to 1 we can 
take the limits tan/3fl —> 1, Z?^ — ^ and find for these 
approximations the upper bound 



IfJ-Rl 



< 



4F|/4 



(40) 



to avoid tachyonic states. Here we have also set Ts — 
as this is numerically always small. For completeness 
we note, that l^s^sl ^^ bounded from above by the re- 
quirement that all couplings stay perturbative up to the 
GUT scale and from below by the requirement of correct 
symmetry breaking as discussed in Sec. |IID[ 

Combining all requirements we find that a light sneu- 
trino, which could be the NLSP, if one of the diagonal Ys 
entries is rather small, < 0.2, and the other two are large, 
^ 0.7. In an abuse of language we call this state a sneu- 
trino even though the corresponding state is dominantly 
a S. However it still can have a sizable i''' admixture 
as exemplified in Fig. M Note that taking Yj"^ small is 
an arbitrary choice ana we could have equally well taken 
one of the two other generations. The smallness of this 



I 

coupling also implies that one of the heavy quasi-Dirac 
neutrinos is significantly lighter than the other two but 
we find that this state is always heavier than the lightest 
sneutrino. Therefore a sneutrino NLSP decays always in- 
visibly into vG. As can be seen in Fig. [o] the next heavier 
state is Xi which turns out to be mainly a xr higgsino. 
If kinematically allowed it will decay to Oifh yielding 



Xi -^ ^i^h 



vGW^*H 



(41) 



giving a final state with an (off-shell) ly-boson, the lep- 
ton of the corresponding generation and missing energy. 
For completeness we note that we find BR{x1) — >• Dij/h 
1 if \Yf 

BR{x°i) 



< 0.07 for the parameters used in Fig. M and 
-)■ f^i z^ ~ 1 for larger values of | Y^^ \ . The latter 
leads to a completely invisible final state and, thus, in 
this part of the parameter space this scenario cannot be 
distinguished from the Xi NLSP case in this model. 



D. Z' phenomenology 

The LHC collaborations ATLAS and CMS have re- 
cently updated the bounds on Mz' from the search for 
dilepton resonances |118l 1119) at ^/s = 8 TeV and an 
integrated luminosity of about 20 fb~^ each. In order to 
apply these bounds to our model we calculate the pro- 
duction cross section of the Z' and the subsequent decay 
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FIG. 10. (j{pp -^ Z' ^f Z+P) as a function of M^/ for the sce- 
narios BLRI (solid line), BLRII (dot-dashed line) and BLRIII 
(dashed line). The red dotted line shows the exclusion limits 
at 95 % CL obtained by ATLAS [TT8] . 



into a pair of leptons as a function of the Z' niasfrl From 
Figure [T0| one can extract the hmits which depend on the 
underlying parameters. In case of BLRIII (dashed line) 
only standard model decay channels for the Z' are open 
leading to a bound of about 2.43 TeV whereas it can 
be reduced to about 2.37 TeV if in addtion decays into 
heavy neutrinos and sneutrinos are allowed as is the case 
of BLRI (full line). This translates into a lower limit on 
vr of about vb. ^ 6.6 TeV. 

As already mentioned, the colored SUSY particles are 
rather heavy in this model in most of the parameter space 
once the constraint on the Higgs mass is imposed imply- 
ing that the discovery of supersymmetry requires either 
a huge statistics and/or a larger c.m.s. energy. However, 
it has been shown that the decays of the Z' open the pos- 
sibility to produce SUSY particles [73 l [77 l fT22Hl26] . As 
has been discussed in [TJ in a CMSSM-like variant of this 
model, the potentially interesting final states from Z' de- 
cays are: v^Vh, II, vv, xtxT ^nd x°x?- However, it turns 
out that in the GMSB variant the required conditions for 
the different channels are harder to realize as this model 
is more constrained. In particular we hardly find charged 
sleptons except for the case that Mz' is above 4 TeV be- 
cause it couples significantly stronger to L-sleptons than 
to _R-sleptons. In Tab. |v]we list the Z' decay modes of 



the parameter points given in Tab. Ill 



The most important non-standard decays of Z' are 
those into heavy neutrinos. Their masses are propor- 
tional to y^Y^ + YgVR implying that the corresponding 
Yukawas should not be too large because otherwise these 
decays are kinematically forbidden. This can clearly be 
seen by combining tables |III| and [V} the smaller Ys the 
larger the corresponding branching is (for fixed Y,^). The 



^ We used CalcHEP 3.4.2 |120| for the cross section calculation. 
The model was implemented using the SUSY Toolbox |121| . 





BLRI BLRII BLRIII 


BLRIV 


BLRV 


BLRVI 


Mz- [TeV] 


2.5 


2.7 


2.4 


4.3 


BR{dd) 


0.45 


0.49 


0.52 


0.52 


0.52 


0.48 


BRiuu) 


0.08 


0.09 


0.10 


0.10 


0.10 


0.09 


BR{ll) 


0.17 


0.18 


0.20 


0.20 


0.20 


0.18 


BR{uv) 


0.15 


0.16 


0.17 


0.17 


0.17 


0.16 


BR{W+W-) 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


BR{vhVh) 


0.12 


0.06 


- 


- 


- 


- 


BR{hiZ) 


- 


- 


0.01 


- 


- 


- 


BR{h2Z) 


- 


- 


- 


0.01 


- 


- 


BR{ir) 


- 


- 


- 


- 


- 


0.02 


BR{i>v) 


0.01 


— 


— 


— 


- 


0.01 


BR[x°X°) 


- 


- 


- 


- 


- 


0.02 


BRixtx2) 


- 


- 


- 


- 


- 


0.02 



TABLE V. Br anch ing ratios of the Z' boson for the parameter 
points of Tab. Ill Only branching ratios larger than 10~^ are 
shown. 



heavy neutrinos decay into Wl, Zv and hv with a branch- 
ing ratio of ~ 0.6, ~ 0.2 and ~ 0.2, respectively |77j. Here 
h denotes the doublet-like Higgs boson. 

Naively one would expect that also sneutrinos should 
show up in such scenarios because as discussed in 
Sec. |III C 3| the smaller Yg the smaller the mass of the 
lightest sneutrino. However, at the same time the S com- 
ponent increases as can be seen in Fig. |9] which reduces 
the coupling to the Z' . For an intermediate range of 
Yg the second lightest sneutrino can be produced in Z' 
decays. It decays dominantly according to 



V2 



yielding 



Vh^dl 



Z' 



VhVhVl -^ IIWW + ^T : 



i>2i>2 



414W^T 



(42) 



(43) 



as final state. The other final states are 212W2ZI^t, 
2l2W2h^T, AZ^Ti'2.Z2h^T and 4/i^t, where h denotes 
again the doublet-like Higgs boson. We note for com- 
pleteness that for some part of the parameter space also 
the decay V2 — >■ i'lh^n is possible offering in principle a 
possiblity to observe /i^^ . As i>2 is relatively heavy, there 
is a kinematical suppression and we find only branching 
ratios of at most 0(0.01) for sneutrinos in the final state. 
The Z' couples to all higgsino states and the corre- 
spond coupling is proportional to 



dxi^iZlfZ^^ 



zi^z]i^) ^ ^{Zl^'Zlf' - Z'^^Zf)) 



(44) 

Here, Z^- is the unitary matrix which diagonalizes the 
neutralino mass matrix. As discussed above, the Xfl-like 
neutralinos can be rather light. However, its admixture 
is such that Z^'^ ~ Z^-'^ ~ ±1/^2 and thus this final 
state has only a tiny branching ratio. For large Mz' the 
decays into to the heavy MSSM charginos/neutralinos 
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containing a sizable higgsino component are kinemati- 
cally allowed with a branching ratio of a few per-cent as 
can be seen in Tab. |V| In this scenario the production 
cross section for the Z' is about 1 fb and, thus, again 
large statistics are needed to observe and study the cor- 
responding final states. 



E. Lepton flavor violation 

So far we have assumed that neutrino mixing is ex- 
plained by the flavor struture of [is which hardly plays 
a role for the phenomenology discussed so far. In this 
case also the rates for flavor violating lepton decays are 
very small and cannot be observed in the near future. 
However, as can be seen from the seesaw approximation 
of the neutrino mass matrix |127j , 



AS 



.,2 ^^ 



^Ys^y.s{Yl)-^Y. ■ 



(45) 



neutrino mixing can also be explained by the flavor struc- 
ture of Yy and Yg ■ As one can always find a basis where 
Yj is diagonal on the expense of having non-diagonal 
y^ and /is we will now consider the other extreme case 
where the complete flavor structure resides in Y^. Non- 
diagonal Yu induces also non-diagonal entries in the soft- 
breaking terms of sleptons. However, as the scale for the 
GMSB boundary is lower than the GUT scale the cor- 
responding effects are significantly smaller compared to 
typical SUGRA scenarios. Therefore, the main effect is 
due to the vertices where the off-diagonal entries of Yy 
enter. To study the corresponding effects in our model 
we parametrize Yy according to [58] : 



y.=] 



/ 

ail--^) _a(l + ^)l, 

2 sin 6*13 1 1 



Ar 



'Atm 



0.4, 



(46) 



using the latest data from the global fit of the PMNS 
matrix jl28j . This fixes Yy up to a global free prefactor 
/ which determines the rate for flavor violating decays 
like /i — > 67 or /i — > 3e. Their branching ratios are con- 
strained by experiment to be smaller than 5.7- 10"^'^ and 
10~^^, respectively |1291ll30j . In addition the /i-e conver- 
sion rate (CR) in gold turns out to be important which 
has to be smaller than 7 • 10" ^^ [HI]. 



In Figure [TT] we show these rates as a function of / 
for the two points BLRI and BLRIII. We observe that in 
these scenarios CR{fj, — e, Au) is the most constraining 
observable followed by BR{fj, — > 3e) and/or BR{ii — > 57). 
As explained in detail in refs. |132j and [13311134] this be- 
haviour can be understood as a non-decoupling effect of 



the Z-boson contribution to CR{fi — e, Au) and BR{fj, — >■ 
3e) which are enhanced by a factor {msusY /Mz)'^ with 
respect to the off-shell photon contributions which is siz- 
able due to the required heavy SUSY spectrum. In all 
cases the sneutrino-chargino loops give the dominant con- 
tributions. We find that CR{ijL~e, Au) gives the strongest 
constraint on the size of / for all points of Tab. |III[ For 
completeness we note, that once the bounds on this ob- 
servable is fulfilled we find that the corresponding LFV 
decays of the r are so suppressed that they are below the 
reach of experiments in the near future. This implies on 
the other hand, that if for example one of the LHC ex- 
periments observes for example r — )■ 3//, then this class 
of models is ruled out. 

A point worth to mention here is a rather strong de- 
pendence of BR(fi — )> 67) on the sign of [Ir which can 
change the rate by an order of magnitude. The reason is 
that the off-diagonal elements for the i'''-S mixing in eq. 



(381 is dominated by the ijrVh contribution as discussed 
in Sec. |IIIC3| We exemplify this behaviour in Fig. [T2| 
where we show BR{fi — >■ ej) as a function of M for both 
signs of fj,ii, f — 0.03 and fixing the other parameters as 
for BLRIV. The black lines give all contributions whereas 
for the light green lines we have taken out the contribu- 
tions containing the higgsino-like chargino. As a conse- 
quence BR(fj, — >■ 67) could be in the reach of an upgraded 
version of the MEG experiments if hr is negative. For 
completeness we note, that this is a specific feature of the 
GMSB model as for example in SUGRA inspired models 
large trilinear couplings Ts could be present dominating 
this mixing. 

Finally we stress again, that the finding of this sec- 
tion depends on the assumption that the complete flavor 
structure needed to explain neutrino data is present in 
Yy. If this structure is at least partially shifted to /is one 
can reduce the predictions for the lepton flavor violating 
observables. It turns out that the size of this reduction 
depends on the SUSY parameters and, thus, we do not 
discuss here in detail. 



IV. CONCLUSION 

We studied in this paper the GMSB variant of a SUSY 
model with an extended gauge sector where the couplings 
unify at the GUT scale. Compared to GMSB with MSSM 
particle content only one can obtain a tree-level mass for 
the light doublet-like Higgs boson above Mz which eases 
the need for large radiative corrections to obtain a Higgs 
mass at 125 GeV. For this reason we find in this model 
a lower bound on the mass of the lighter stop of about 
two TeV which is about a factor two smaller than in 
usual GMSB models. Nevertheless the SUSY particles 
of the strongly interacting sector are rather heavy and, 
thus, the existing bounds on squarks and gluino are satis- 
fied automatically. However, this implies that one needs 
high luminosity at the LHC with y/s = 14 TeV to study 
this sector and the resulting cascade decays in detail. 
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FIG. 11. Flavor violating observables as a function of / as defined in eq. (46 1: BR{iJ, — >■ 67) (solid line), BR{fi — >■ 3e) 
(dotted line) and CR{fi — > e) in Au (dashed line) for the points BLRI (left) and BLRIII (right) defined in Tab. Ill and using 
Y^ to explain the neutrino data. The upper bounds {BR{fj, -i> e-y) < 5.7 • 10"" [H^ . BR{^j. -> 3e) < 1.0 ■ 10"^^ 
CR{fj. — e, Au) < 7.0 • 10~^^ |131| ) are shown as a red horizontal line, respectively. 




1 X 10'" 2 X 10'" 5 X 10'" 1 X 10" 2 X 10" 



M [GeV] 



FIG. 12. Branching ratios of the LFV decay BR{jj. — >■ 67) as 
a function of the messenger scale M for sign {jj.r) — + (black 
solid line) and sign {fin) — — (black dashed line) fixing the 
other parameters as in BLRIV and / = 0.03. The respective 
straight green (light) lines show the branching ratio excluding 
the contribution of charged higgsinos and sneutrinos. 



what harder to achieve than in usual GMSB models. We 
have worked out main features of the corresponding sce- 
narios paying also particular attention to the possiblity 
that the new Higgs boson, which can be rather light, can 
show up in the SUSY cascade decays. We have argued 
that Z' decays can serve as SUSY discovery even in this 
rather restricted model. 

Last but not least we have discussed which lepton fla- 
vor violating observables can be observed in this class of 
GMSB models. It turns out that /z — >■ 3e and fj, — e- 
conversion are usually more constraining than fi — ;> ej. 
The necessary requirement for a possible observation is 
that there are sizable off-diagonal entries in the neutrino 
Yukawa coupling. It turns out that the rates for the cor- 
responding r-decays is well below the sensitivity once the 
contraints from the niuon-sector are taken into account. 
Note, however, that the rates for all observables get tiny 
if neutrino mixings are explained via the flavor structure 
of ij.g instead of the flavor structure of Y^, . 
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Appendix A: Mass matrices in the U{1)r x U{1)b-l 
basis 

Here we give for completeness the mass matrices that 
were shown in the text for the original basis of SU{3)c x 

SUi2)LXU{l)RXUil)B^L- 



1. Higgs mass matrix 



In the basis (ad, o'u, (Jr, or), the Higgs mass matrix is given by 



m%, 



( \9iL^'^c\^m\s^ 



S2/3 



\g'^vvRCfiCfi^ -\g'^vvRSpCfjj^ 

-Yg'^VVRCpSp^ \g'^VVRSi3Si3„ 



\g'^vvRCfjCii^ 
-\9^vvRSfiCf,^ 

1 x2 „,2 „2 



-^VVRCpSp^ 
Ig'^VVRSftSp^ 

*2^H //l™2 I ~2 



I~9rR^rcL + ^Ar<4n -^(4"^i« + ~9rR^r) 

- ^ (4™^ + 9JiRvl) WrrvIsI^ + m\^cl^ j 



16 



(Al) 



where g^ = gR{gR - gSLR) + gRBL^gRBL - gBL), CiIl = 9l+9r + 9rbL and g|j^ = {gBLR - gR^ + {gBL - gRBbf- 



2. Neutralino mass matrix 



The neutralino mass matrix in the basis (A^-l, A^/a, /ij], /i", A;?, Xi?,Xi^) reads 



M^o ^ 



V- 



Mb-l 


gRBLVd 
2 
gRBLVu 
2 

Mblr 

(gBL-gBBL)VxR 
2 
igBL-gRBL)VxR 

2 





ghVd 

2 

gLVu 

2 







gBBLVd gRBLVu 

2 2 

2 2 





gnVd 
2 









2 









2 

gRVu 
2 

M, 



{gR-gBLR)vxR 

2 

JgR-gBLRjVxR 

2 



gBL-gRBL)VxR {gBL-gRBL)VxR \ 

2 2 1 







(gR-gBLR)'"xR {gR-gBLR)VxR 







2 





(A2) 



3. Slepton mass matrix 

In the basis (e^, cr), the slepton mass matrix is given by 



m? 



™fL 



72 



(T]cos^-;uy^^sin^)' 



%{T,cosp~^l*Y,sm|i) 



.V2 



^L 



(A3) 



with 



™[L 



=™i + 2^^ *^"^^ PY^^Ye + o ( ^ (.9l - 9R9BLR - gBLgRBL)v^ COS 2;3 

^ (sIl + 9%LR - 9R9BLR - 9bl9rbl)vr cos 2/3ij ) 1 , 



"^Ffl 



L "™e'= + 2^^ '^"'^^ PYeY^ ^ 8 ( ^ {9r{9BLR + 5fl) + 9RBL{gBL + gRBL)v'^ COS 2/3 



(sIl - .91 + g'BLR - gRBL)vR COS 2/3^ ) 1 . 



(A4) 



Sneutrino mass matrix 



The sneutrino mass matrix in the basis {v, i'^, S) reads 



/ 



M? = 



ml, 



L V2 

vr 



UTlsp-^^Y^p) 



V2 
1 



Vfl 



V ^vvRYlY^spSfi^ ^{Tlsfs^-^RY^Cfs^) 



^vvrYJYss^spj, \ 

(A's+Ms)(ms+A's 



2 



^^i>S 



(A5) 



with 
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m 



vL 



'■OR 






g \y'^C2i}{gi + gRgBLR. + 9bl9rbl) - «flC2/3„(5lL + 5sLfl - 9R9BLR " 9bl9rbl)) 1 : 
5 + y^^'^^'^^ + 8 («^C2/3((g/? - 9blr)9r + 9rbl{9rbl - 9bl)) + 



^fl„2 



^fe^S^5 



^^^^2/3^ ((5fl - 9BLr)^ + {9BL - 9RBLf) ] 1 • 



(A6) 
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